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Q-i ABSTRACT 

q . Short gamma-ray bursts (SGRBs) observed by Swift are potentially revealing the first 

$H 1 insight into cataclysmic compact object mergers. To ultimately acquire a fundamental 

I understanding of these events requires pan-spectral observations and knowledge of 

their spatial distribution to differentiate between proposed progenitor populations. 
Presently (late 2011) there are only some 30% of SGRBs with reasonably firm redshifts, 
and this sample is highly biased by the limited sensitivity of Swift to detect SGRBs. We 
account for the dominant biases to calculate a realistic SGRB rate density out to z « 
\ 0.5 using the Swift sample of peak fluxes, redshifts, and those SGRBs with a beaming 

[ — . angle constraint from X-ray/optical observations. We find an SGRB lower rate density 

of 7.1^3 2 Gpc _3 yr~ 1 (assuming isotropic emission), and a beaming corrected upper 
limit of 1200±j^ Gpc 3 yr 1 . Assuming a significant fraction of binary neutron star 
CN \ mergers produce SGRBs, we calculate lower and upper detection rate limits of (1 — 

C> 200) yr -1 by an ALIGO and Virgo coincidence search. Our detection rate is similar 

to the lower and realistic rates inferred from extrapolations using Galactic pulsar 
observations and population synthesis. 



Key words: stars - gamma-ray burst: individual - gravitational waves - techniques: 
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1 INTRODUCTION 

It is generally accepted that 'long' (T90 > 2 s) gamma- 
ray bursts (GRBs) are l inked to the core collapse of 
massive stars (collapsars') |Wooslevl 1 19931 : iPaczvnskl 1 19981 : 
iMacFadven fc WooslevH i999'). For several cases, the GRBs 
are firmly associated with Type Ib /c supernovae (e.g. 
iHiorth et al]|2003l : IsTanek et al.l l2003h . This strongly sug- 
gests that they are linked to the end point of massive stellar 
evolution. In contrast 'short' (T90 < 2 s) gamma-ray bursts 
(SGRBs) have a less certain origin. 

The first breakthrough to understanding the nature of 
SGRBs was made in 2005 after the launch of the NASA 
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1 T90 is the duration in which the cumulative counts are from 5% 
to 95% above background. 



Swift satellite (|Gehrels et all 12004 ). Prompt localizations 
and deep afterglow searches yielded the first redshifts and in- 
vestigations of their progenitor environments. By late 2011, 
about three dozen SGRBs had been localized by Swift. 
Among them, about 50% have optical detections and about 
one third have redshift determinations based on host galaxy 
spectroscopy. 

Binary neutron star mergers (NS-NS) or neutron star- 
black hole (NS-BH) mergers are favoured as the pro- 
genitors for SGRBs, based on the association of some 
SGRBs with an older stellar population, as compared to 
long GRBs. Evidence for the origin of SGRBs in the fi- 
nal merger stage comes from the host galaxy types (e.g. 
iLee Ramirez-Ruiz fc Granodl2005lZheng fc Ramirez- Ruiz I 
120071 ). Kicks imparted to NSs at birth will produce velocities 
of several hundred km s -1 , implying that binary inspiraling 
systems may occur far from their site of origin. Fong, Berger 
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& Fox (2009) using Hubble Space Telescope observations to 
measure SGRB-galaxy offsets, find the offset distribution 
compares favourably with the predicted distribution for NS- 
NS binaries. 

Despite this progress, unambiguous identification of 
burst types continues to be a problem. The various am- 
biguities have motivated several authors to re-define dif- 
ferent classes of GRBs via a number of properties includ- 
ing: spectral features, associated supernova, stellar popula- 
tion, host galaxy, location in the host galaxy and progen- 
itor type (Zhang et al l 120071 ) . Some authors also suggest a 
third interm ediate population of bursts based on the existing 
Tgo scheme (|Horvath et al.ll2008l) with a lower than average 
peak-flux distribution jveres et alj|2010l). 

Furthermore, work by iBromberg et al.l l|201ll ) shows 
that categorizing bursts as SGRBs based on Tgo is satellite 
dependent. They find that 2 s > Tgo is statistically reason- 
able for BATSE bursts, but not accurate for Swift bursts. 
To quantify this uncertainty, for a 0.7-s burst they show 
that there is an equal probability of its being either long or 
short. We account for this uncertainty, and its effect on the 
intrinsic SGRB rate, by using their probabilities to scale the 
calculated SGRB rate. 

Around 20% of the SGRBs detected by Swift have been 
followed by an ex tended emission lasting up to 100 s (here - 
after SGRB-EE) dNorris fc BonnelfeOOrirPerlev et al.ll2008h 
leading to suggestio ns that different progenitor typ es pro- 
duce these bursts (jNorris. Gehrels fc Scargldlioilh . Can- 
didate systems include the birth of a rapidly rotating 
proto-magnetar produced via NS-NS merger or accretion- 
induced collapse of a white dwarf (jMetzger et al.l 120081 ; 
iBucciantini et al] l201lh ■ In the following, we have consid- 
ered separately the normal SGRBs and the SGRB-EE to 
take this uncertainty in the progenitor nature into account. 

SGRB observations could potentially play a more im- 
portant role in constraining compact object merger rates on 
two fronts as more data become available. Firstly, as a di- 
rect consistency check of the binary merger rates inferred 
from binary pulsars and population synthesis. Secondly, as 
a constraint on the progenitor scenario and the engine that 
drives the emissions. Both are critical issues for solving the 
puzzle of the origin of SGRBs and their links to compact 
object mergers. 

There are two popular but different methods employed 
for estimating NS-NS merger rates. The first uses extrap- 
olation from the observed sample of NS binaries infered 
from pulsar observa tions (Naravan. Piran fc SheniH Il99ll; 
iKalogera et~aH 12006): the second uses population synthe- 
sis simulations ( Belczvnski et all [2007). where several un- 
known model parameters are constra ined by obser v ations 
and others are assumed from theory. lAbadie et ail <|2010l ) 
(hereafter AblO) assessed the rates from these techniques, 
defining lower and realistic rate densities of 10 and 1000 
Gpc _3 yr _1 respectively. These rates are uncertain by 1 — 2 
orders of magnitude because of the small number of ob- 
served Galactic binary pulsars, and from poor constraints 
in population-synthesis models. 

Within 5 years, co-ordinated gamma-ray, X-ray, optical 
and gravitational-wave observations may allow the strong 
gravity regime of the central engine of compact object merg- 
ers to be probed. Such 'multi-messenger' observations pro- 
vide the opportunity to probe these events across a vast 



energy spectrum, and to constrain the progenitor popu- 
lations of SGRBs. Furthermore, co-ordinated optical and 
gravitational-wave searches may play an important role in 
confirming the first direct gravitational-wave o bservations of 
compact object mergers (jCoward et a.1.1 [201 lh ■ It is becom- 
ing increasingly important to constrain the rate of compact 
object mergers and their proposed optical counterparts in 
the context of up-coming gravitational-wave searches. 

The scarcity of SGRB observations and the poorly un- 
derstood biases in the current data have made it challeng- 
ing to calculate the SGRB rate density with meaningful 
uncertainties. Previous work using GRB flux-limited sam- 
ples either assumes, or derives models for, the SGRB lu- 
minosity function and rate evolution of the sources; e.g. 
iGuetta fc Piranl i|2006h find a SGRB rate density of 8 - 30 
Gpc~ 3 yr~\ assuming isotropic emmission. Another study 
(|Dietj201ll ). calculates a higher NS-NS merger rate of about 
7800 Gpc^V- 1 , but assuming the emissions are beamed 
with small opening angles. 

In this Letter, we calculate a beaming-corrected SGRB 
rate density using the Swift sample of SGRB peak fluxes, 
redshifts and inferred beaming angles from X-ray observa- 
tions. For our GRB selection criteria we use the Jochen 
Greiner catalogue of localized GRBs (see Table 1) and se- 
lect bursts indicated as short that have reliable redshifts up 
to 2011 August. From this selection of 11 bursts, we omit 
the burst GRB 090426, linked to the death of a massive 
star (as discussed in section 2). We note also that although 
GRB 100816A has shown certain characteristics not consis- 
tent with a SGRB classification it cannot be clearly ruled 
out; additionally, its redshift of 0.8 will not influence our 
results significantly. 

We avoid using a SGRB luminosity function, models for 
progenitor rate evolution, and a beaming angle distribution, 
all of which have large uncertainties. Instead, we focus on 
observed and measured parameters that take into account 
selection effects that modify Swift's detection sensitivity to 
SGRBs. Finally, we use our SGRB rate density estimates 
to infer a detection rate of binary NS mergers by Advanced 
LIGO (ALIGO) and Virgo interferomters. Despite the poor 
statistics, this approach gives meaningful results and can 
be followed up when a larger sample of SGRB observations 
becomes available. 



2 SGRB BEAMING CONSTRAINTS 

The currently favoured emission model for SGRBs 

is a compact object merger triggering an ex- 

cau si ng a burst of collimated 7-rays 



plosion 



7-rays 

ilEichler et all Il989l; iNarayan. Paczvnski fc Piranl Il992l ; 
iLee Ramirez-Ruiz fc GranotJ 20051 ) powered by accretion 



onto the newly formed compact object. The ultra-relavistic 
outflow is eventually decelerated by interaction with the 
interstellar medium to produce a fading X-ray and op- 
tical afterglow. After the jet Lorentz factor F decreases 
to r ~ 0T , where 6*j is the jet opening half angle, the 
afterglow becomes observable from viewing angles greater 
than Oy 

Only two SGR Bs have identifiable jet breaks (see 
ICoward et aT]|201ll . and references therein): GRB 050709 
(6j ~ 14° - detected by the HETE satellite) and GRB 
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051221 (0j ~ 7°). These beaming angles are not necessarily 
agreed upon, mainly because the afterglow jet break times 
obtained from X-ray or optical light curves are notoriously 
difficult to identify at late times, and the emissions are gen- 
erally fainter in X-rays and optical compared to those of 
long bursts. Furth ermore, the inferred beaming angles are 
model dependent (|Sari. Piran Sz Hafpern 1 1 1998T ) and sensi- 
tive to the circumburst environment. 

For the following rate density calculations, we employ 
the standard fireball model (e.g. Sari et al. 1999) to estimate 
the jet opening half angle using the minimum jet break times 
from X-ray observations: 



9j = 0.161^/(1 + z)] 3/8 (nr 1 /E iao ) 1/s rad, 



(1) 



where tj, the jet break time, is the time in days when the 
light curve decay index changes from —1 to —2. The param- 
eter n is the circum-burst density, in cm -3 , r\ is the fraction 
of GRB energy observed in the prompt emission, and -Ei ao is 
the prompt isotropic equivalent energy in 10 52 erg. We as- 
sume an average of n = 1 cm" 3 and r) — 0.1 for our SGRB 
sample. 

The difficulty of measuring jet breaks leads to a bias in 
obtaining jet angles favouring those bursts that have very 
bright X-ray/optical afterglows and/or relatively short jet 
break times. Hence, the true average of the SGRB jet an- 
gle distribution will be greater than that based on the very 
small sample used in this work. Nonetheless, given the lack 
of knowledge of the SGRB angle distribution, both theoret- 
ically and observationally, we choose to employ the small 
number of inferred SGRB beaming angles and lower limits 
from X-ray /optical observations. 

To extend the very small measured beaming angle sam- 
ple, we constrain the minimum beaming angles for those 
SGRBs that have X-ray afterglow light curves available up 
to at least 1 day from the trigger. The following five SGRBs: 
GRB 100816A, GRB 071227A, GRB 070714B, GRB 061006 
and GRB 050724 satisfy this criterion. We set a lower limit 
on the jet opening angle from the time when the Swift XRT 
monitoring stopped. We compute Ei so , from the literature, 
if available, or from the quoted fluences in the GCN circu- 
lars by estimating the corresponding Ei BO in the 1-10 4 keV 
rest-frame energy range. For the GRB energy spectru m, we 
use t he broken power law fit, termed the Band model |Bandl 
2006), or cut-off power law with spectral parameters where 
known. If the spectral parameters are unkown we simply 
rescale the fluence by 47rdf,/(l + z), where <Il is the lumi- 
nosity distance. 



3 EMPIRICAL SGRB RATE MODEL 

A Poisson GRB rate can be estimated from small numbe r 
statistics using V max jPiran I Il992l ; ICohen fc Piranl 1 19951 ). 
where V^ax is the maximum volume within which a SGRB 
with observed peak flux, F p , could be detected for a given 
satellite detector with detection sensitivity -FLim- While I/ ma x 
has been used previously (e.g. iGuetta fc Delia Valle] 2007 ) 
to es timate rates of low- luminosity GRBs (see also ICowardl 
2005), we extend the rate estimate method to account for 
biases in the Swift GRB redshift and peak flux sample. 

Firstly, the low energy detection bandwidth of Swift 



SGRB 


T90 




z 


peak flux 


rcf 




(s) 


(deg) 




(ph s _1 cm -2 ) 




101219A* 


0.6 


- 


0.718 


4.1 


(1) 


100816A* 


2.9 


> 12T 


0.803 


10.9 




100117A* 


0.3 




0.92 


2.9 


(3) 


090510A* 


0.3 




0.903 


9.7 


(4) 


080905A 


1.0 




0.122 


6.0 


(5) 


070724A 


0.1 


> lit 


0.457 


2.0 


(6) 


061217A 


0.2 




0.827 


2.4 


(7) 


060502B 


0.1 




0.287 


8.5 


(8) 


051221A* 


1.4 


7 a 


0.547 


12.0 


(9) 


050509B 


0.73 




0.225 


3.7 


(10) 


SGRB-EE 


071227A 


1.8 


> 16+ 


0.383 


4.4 


(11) 


070714B 


64 


>6t 


0.923 


10.0 


(12) 


061210 


85 


> 12+ 


0.409 


62.9 


(13) 


061006 


129 


>6t 


0.437 


15.8 


(14) 


060614 


109 


12 b 


0.125 


25.0 


(15) 


050724 


96 


>25t 


0.258 


12.8 


(16),(17) 



Table 1. SGRB peak fluxes, T 90 and red- 
shifts taken from the Swift online catalogue and 
http://www.mpe.mpg.de/~jcg/grbgcn.html used to calcu- 
late Poisson rates. We use the 20-ms peak photon fluxes from 
the BAT2 catalogue where possible - those marked by * are 1-s 
peak photon fluxes, t We set a lower limit on the jet opening 
angle from the time when the Swift XRT monitoring stopped. 
a The detection of a jet break was observed for GRB 051221A 
in X-rays (Soderberg et al. 2006) at about 5 days post-burst, or 
6 ~ 7°. b jet break identified at 1.4 days post-burst, circum-burst 
density n = 10 cm -3 , and E- lso = 2.5 X 10 51 ergs. 
Refere nces -(l) GCN 11518, (2) GNC 11123, (3) iFong et all 
||2011|). f4) lAckermann et all ||2010|). (5)]Rowlinson et alj 120101), 
(6) IBerger et alj ||2009|). (7) IBerger et all ll2007l), (8)lBloom et all 
||2007|) . (9) ISoderberg et all j2006h . (10) iGehrels et all <2005l) . 
The SGRB s below the line are classified as SGR B-EE - (11) 
ICaito et all d201Cl). (12)lGraham et all j2009l). (13) IBerger et all 
||2007|). (14) iBerger et all ll2007h. (15) bella Valid ll2006h , (16) 
iBerger et all j2005l ). (17) iGruppe etall l|2006l ). 



(15-150 keV) in comparison with BATSE's 50-300 keV re- 
sults in a bias against SGRBs which typically have harder 
emissions. Secondly, the Swift detection threshold is not sim- 
ply defined by the detector sensitivity, but by a complex 
triggering algorithm. Both these effects manifest as a bias 
against Swift detecting SGRBs; i.e. a smaller proportion of 
bursts has been detected by Swift (~10%). The latter effect 
results from the detection process employed by BAT, the 
Swift coded-aperture mask 7-ray detector. In addition to 
requiring an increased photon count rate above background 
(the sole triggering criterion used for BATSE) , BAT employs 
a second stage in whi ch an i mage is formed by accumulating 
counts for up to 26s lBandl (|2006l ). 

We attempt to crudely correct for these biases by using 
the observed SGRB rate from BATSE as a rate calibration 
for the Swift SGRB rate. Because BATSE operated at dif- 
ferent energy thresholds and trigger sensitivities, for consis- 
tency we take all BATSE SGRBs with 64-ms peak flux when 
the trigger threshold was set to 5.5 a in the 50-300 keV en- 
ergy range (total live operation time of 3.5 years). This yields 
32 S GRBs sr" 1 y r" 1 , assuming an effective BATSE FoV of 
n sr iBandl (|20031 ). The ratio of the rate of BATSE to Swift 
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SGRBs, Rb/s = 6.7, is used to calibrate the observed Swift 
SGRB rate to estimate an intrinsic SGRB rate. 

Because SGRBs occur over a short duration, it is more 
difficult (compared to long bursts) to produce a significant 
signal above background. Hence instead of using the theoret- 
ical BAT sensitivity of FLim = 0.4 phs" 1 cm -2 , we employ a 
flux limit of 1.5 phs -1 cm -2 , using the smallest 20-ms peak 
flux from the SGRB data. 

Another detector-dependent selection effect arises from 
the limited energy bandpass of Swift, i.e. (15-150) keV. 
What is measured by Swift is not the bolometric peak flux, 
but a detector response and source-spectrum dependent 
peak flux. We define the b olometric isotropic-e quivalent lu- 
minosity correction factor (|lmerito et al.ll2009f ): 



Cdet(ei, e 2 ) 



EN(E)dE 
£*EN(E)dE 



(2) 



where [E\ = l,i?2 = 10000] keV spans the bolometric 
gamma-ray spectrum, and [ei , 62] is the sensitivity band of 
Swift, i.e. 15 to 150 keV. For the s ource energ y spectrum, 
N(E), we employ the Band function (Band 2006), with spec- 
tral indices a = — 1, ft = —2.3 and a rest-frame power-law 
break energy of 511 keV. 

In addition, for high redshift sources, a higher energy 
component of the source spectrum is redshifted into the sen- 
sitivity band of the detector. The following factor, k(z), ac- 
counts for the downshift of 7-ray energy from the burst to 
the observer's reference frame: 



k(z) 



g EN(E)dE 



(3) 



Taking into account the sensitivity reduction and k- 
correction, the SGRB all-sky rate can be inferred from the 
flux-limited SGRB sample in TableQ] We calculate the max- 
imum distance d max , with corresponding redshift z max , that 
a burst at luminosity distance dh could be detected given 
Swift's sensitivity FLim: 



fin 



F p k{z) 



-FLimfc(zLim) 



d L {z) 



(4) 



where F v is the observed peak flux and fc(zLim) is the k- 
correction for a burst at the maximum limiting distance of 
detection for the sample. Given the large uncertainty in the 
SGRB luminosity function, we use the largest redshift from 
our sample, ZLim = 0.92. We note that the bolometric correc- 
tion, Cd e t(ei, e2), cancels out in the above equation, because 
it is independent of redshift. 

The corresponding maximum SGRB detection volume 
for each burst is defined as 



V m 



dV 
dz 



dz , 



(5) 



where the volume element factor, dV/dz, and luminosity dis- 
tance, dL(z), are calculated using a flat- A cosmology with 
H = 71 km s" 1 Mpc -1 , U M = 0.3 and fi A = 0.7. This 
estimate (Eq. 5) is valid as long as the rate of SGRBs does 
not change significantly over the range (0, z max ). Since the 
detection range of SGRBs is quite small and typically z < 1 
this approximation is valid. 

To calculate the intrinsic rate of SGRBs requires ac- 



counting for the beaming angle, 9j , of the jetted burst. Equa- 
tion ((6j) expresses the beaming factor used to correct for the 
unobserved SGRBs that are not detected because the jet is 
misaligned with the detector: 



B(e j ) = [l-cos(e j )]- 1 



(6) 



To account for the fact that only a fraction of observed 
SGRB have measured redshifts, we scale the rate density 
by the ratio of Swift bursts with redshift to those without 
redshifts, F r . We use F r « 10/41 and 6/12 for SGRB and 
SGRB-EE respectively. The time span encompassing all ob- 
servations, T ~ 6 yrs, is defined by the start of Swift obser- 
vations to the time of the most recent SGRB in the sam- 
ple and we account for the fractional sky coverage of Swift, 
Q « 0.17. To account for Swift's reduced sensitivity for de- 
tecting SGRBs relative to BATSE, we use the ratio of the 
BATSE to Swift SGRB detection rate, which we approxi- 
mate as Rb/s — 6.7. This converts the Swift SGRB rate to 
an intrinsic SGRB rate. We point out that this correction 
applies only to SGRB, and not SGRB-EE, because Swift is 
more sensitive to these longer emission bursts. Finally, we 
scale the rate density by Pi(Tgo), the probability that the 
burst is an SGRB based on its Tgo. Combining all detection 
parameters yields the Poisson SGRB rate density of a single 
(ith) burst, and the total rate density for n bursts: 



Rsgi 



E 



1 1 1 



FrTD. 



R B/ sBi(0 s )Pi(Tu 



(7) 



We calculate lower rate estimates assuming no beaming. 
Upper rates are estimated using the observed beaming angle 
constraints for each burst, and if a beaming angle estimate 
is unavailable, we use the smallest inferred jet angle from 
the data, 9j ~ 6°. Given that the inferred average beaming 
angle of the jet is biased towards SGRBs with short jet- 
break times (see Sj2}, our upper rate limit estimate is at the 
limit of plausibility. 

We use the above SGRB rate density to infer a detec- 
tion rate of binary NS mergers by advanced gravitational- 
wave interferometers. For ALIGO and Virgo interferometer 
sensitivities, the horizon distance, Dh (all sky locations and 
orientation averaged over) for optimal detection of a NS-NS 
merger in a coincidence search is about 340 Mpc (cosmo- 
logical redshift not included). For a direct comparison with 
AblO, the detection rate is computed considering the noise 
power spectral spectral density of a single interferometer, 
with Dh = 197 Mpc. Given the uncertainty in the beaming 
angle distribution, we define an optimal detection rate as a 
function of 9j using the SGRB lower rate estimate, i.e. from 
9j = 90°, scaled by B(9j) and the Euclidean volume: 



—Df l Rj jOW B(9j) 



(8) 



4 RESULTS 

We have focussed on the sensitivity of Swift for detecting 
SGRBs, and highlight several important issues. The com- 
plex triggering algorithm of the BAT is biasing detection 
against SGRBs, and this should be considered when at- 
tempting to estimate an intrinsic rate density of SGRBs. 
The rate estimates using V ma ,x are sensitive to the flux limit 
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SGRB Pi(Tgo) lower rate upper rate 

Gpc -3 yr -1 Gpc -3 yr _] 




10 20 30 40 50 60 70 80 90 
SGRB beaming angle (degrees) 



101219A 


0.79 


0.033 


6* 


100816A 


0.51 


0.01 


0.48 


100117A 


0.89 


0.033 


6* 


090510 


0.89 


0.017 


3.1* 


080905A 


0.72 


4.1 


750* 


070724A 


0.84 


0.65 


35 


061217 


0.89 


0.17 


31* 


060502B 


0.72 


0.35 


64* 


051221A 


0.72 


0.022 


3 


050509B 


0.79 


1.7 


310* 


Total rate 




7 1+ 4 ' 9 
'• 1 -3.2 


1200i^ 


SGRB- EE 


071227 




0.038 


0.97 


070714B 




0.0036 


0.66 


061210 




0.0034 


0.16 


061006 




0.0088 


1.6 


060614 




0.073 


3.3 


050724 




0.03 


0.32 


Total rate 




16+ 015 
u ' lu -0.088 


7 i+6.9 
f.l_ 4 



Table 2. The beaming-corrected SGRB rate densities with Pois- 
son uncertainties using the observed constraints on £>j , and scaled 
by the probability Pi (Tgo) that a given burst is a non-collapsar 
(SGRB) based on its Tgo. Lower rate estimates assume isotropic 
emission, and upper rates use the observed beaming angle con- 
straints shown in Table 1, or the smallest observed beaming angle 
in the sample*, 0j ^ 6°. We calculate the SGRB-EE rate density 
separately, noting that Pi (Tgo ) are not shown because Tgo is gen- 
erally much larger than 2 s, as shown in Table 1. 



of the detector, but are mostly invariant to the bolometric 
and k-correction because these effects cancel out. 

With the above caveats, Table 2 shows the rate densities 
of SGRB and SGR-EE using equation applied to each 
burst, with the observed constraints on 6 y The total lower 
and upper rate densities for SGRBs are 7.1^3 2 (assum- 
ing isotropic emission) and I2OOI550 Gpc^yr^ 1 (assuming 
beaming-see Table 2) respectively, w here the errors are the 
90% Poisson statistics <|Gehrelslll986h . The SGRB rate den- 
sity is dominated by GRB 080905A, which was relatively 
faint and nearby, with a z m ax ~ 0.5. We find a correspond- 
ing rate density for SGRB-EEs of 0.16lg;J| 8 and 7.l±l* 
Gpc _3 yr _1 respectively. It is interesting that the SGRB-EE 
rate density from this work is similar to the expected BH-NS 
merger rates from population synthesis models (AblO). 

Assuming a signicant fraction of binary NS mergers pro- 
duce SGRBs, it is interesting to compare these rate esti- 
mates to the NS-NS merger rate density estimated from 
other work. Our lower rate density that assumes isotropic 
e mission is comparable t o the isotropic emission estimates 
of iGuetta fc Piranl (|2006h i.e 8-30 Gpc" 3 yr _1 . Our upper 
limit is comparable to their estimates of beamed emission, 
240 - 1500 Gpc" 3 yr _1 . Compared to AblO, our limits of 
(7 — 1200) Gpc~ 3 yr _1 are similar to their low-realistic rate 
densities (10 - 1000) Gpc~ 3 yr~\ 

Finally, we estimate a plausible detection rate by 
ALIGO/ Virgo using the SGRB rate density limits and equa- 



Figure 1. The lower curve plots equation (0, the detection rate 
of binary NS mergers by a single ALIGO interferometer (Dh = 
197 Mpc) as a function of SGRB beaming angle, using -Rlow =7.1 
Gpc- 3 yr- 1 (see Table 2). The upper curve assumes an optimal 
coincidence search with ALIGO and Virgo interferometers with 
D h = 341 Mpc. 

tion ((8|. The corresponding lower and upper detection rates 
using the SGRB rates from Table 2 are (0.2 — 40) yr" 1 for a 
single interferometer with Dh = 197 Mpc. AblO give a detec- 
tion rate for ALIGO of (0.4, 40, 400) yr -1 for low, realistic 
and high detection rates respectively. For an optimal coinci- 
dence search at ALIGO/Virgo sensitivities, with Dh = 341 
Mpc, the detection rate increases to (1 — 200) yr -1 . Figure 
1 plots the detection rate for a single ALIGO interferometer 
(Dh = 197 Mpc), and an optimal coincidence search with 
ALIGO/Virgo (D h = 341 Mpc), as a function of SGRB 
beaming angle using equation (JHJ. We use the SGRB rate 
density 7.1 Gpc -3 yr -1 (see Table 2) that assumes isotropic 
emission. 

Our detection rate is similar to the lower and realistic 
rates based on population synthesis simulations using Galac- 
tic pulsar observations. The rates are also compatible with 
those of Guetta & Piran (2006), who employ a luminosity 
function and SGRB rate evolution. 



5 DISCUSSION 

It has been suggested that the observed SGRB population is 
presently not a useful constraint on rate estimates for com- 
pact binary mergers (AblO). The justification of this argu- 
ment stems from the fraction (if any) of binary NS mergers 
that will produce an SGRB. Another issue is that SGRB ob- 
servations by Swift and subsequent follow-up in optical are 
affected by selection effects and uncertainty in the beaming 
angles. If SGRB/SGRB-EE have an intrinsic average beam- 
ing angle of ~ 50°, then the jet-break times may not occur 
until weeks after the burst, when the afterglow is too faint 
in the optical. 

The SGRB rate density depends on the spatial- 
distribution of the sources, but because redshift measure- 
ment depends on a bright afterglow, the spatial distri- 
bution is bias e d to the brighter b u rsts that are nearb y 
|Coward et all |200S| ; ICowardl |2009| ; llmerito et all |2009h . 
Hence, for the majority of bursts with a faint afterglow, red- 
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shift measurement is more difficult. We attempt to crudely 
correct for this bias by boosting the calculated rate by the 
fraction of bursts without measured redshift to those with 
redshift, assuming the missing redshifts follow the same dis- 
tribution as the observed redshifts. This assumption is rea- 
sonable (in this work) because the rate density is dominated 
by those faint bursts in the smallest V ma x. 

Hence, there are two main sources of selection effects 
and biases. Firstly, there are the satellite detectors, exempli- 
fied by the significant difference between BATSE and Swift 
SGRB detection efficiencies. Secondly, there is the problem 
of obtaining a redshift either from the host galaxy or the 
afterglow itself (as highlighted above). We have shown that 
these detection biases should be considered when attempt- 
ing to use the current (and future) SGRB detections for con- 
straining rate evolution, rate densities and linking SGRBs 
to binary NS mergers. 

The binary NS gravitational-wave detection rate esti- 
mates are based on calculating an intrinsic SGRB rate den- 
sity using Swift localized bursts, taking into account domi- 
nant selection effects. This approach, based on observational 
data is very different from that based on Galactic binary pul- 
sar observations and modelled population synthesis. In the 
latter, AblO use the observed Galactic binary pulsar popu- 
lation and extrapolate a NS merger rate density out to the 
ALIGO and Virgo detection horizon. Conversely, our ap- 
proach avoids this extrapolation because it is essentially an 
observed rate extending well beyond the average sensitivity 
distances of the upcoming gravitational-wave searches for 
compact binaries (about 300 Mpc or z = 0.07 for ALIGO 
and Virgo interferometers) and, moreover, a significant frac- 
tion of SGRBs are observed in association with evolved stel- 
lar populations. 

In conclusion, the upcoming gravitational-wave detec- 
tion era will be fundamental for resolving the SGRB-binary 
NS merger connection, since an unequivocal association be- 
tween SGRBs and binary NS mergers will only be possible 
via coincident gravitational-wave and electromagnetic obser- 
vations. Ultimately, a comparison beteween the SGRB rate 
density and the gravitational-wave detection rate will help 
constrain the fraction of binary NS mergers that give rise to 
SGRBs and the SGRB beaming angle distribution. 
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